A novel approach to Isoscaling: the role of the order parameter 
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Abstract 

Isoscaling is derived within a recently proposed modified Fisher model where the free energy 
near the critical point is described by the Landau 0(m G ) theory. In this model m = N ^ z is the 
order parameter, a consequence of (one of) the symmetries of the nuclear Hamiltonian. Within this 

o 
o 



framework we show that isoscaling depends mainly on this order parameter through the 'external 
(conjugate) field' H. The external field is just given by the difference in chemical potentials of 
the neutrons and protons of the two sources. To distinguish from previously employed isoscaling 
relationships, this approach is dubbed: m — scaling. We discuss the relationship between this 
framework and the standard isoscaling formalism and point out some substantial differences in 
interpretation of experimental results which might result. These should be investigated further 
both theoretically and experimentally. 

PACS numbers: 21.65.Ef, 24.10.-i, 24.10.Pa,25.70.Gh, 25.70.Pq 
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INTRODUCTION 



In near Fermi energy heavy ion collisions fragments are copiously produced. The mass 
distributions of these fragments often exhibit a power law behavior 

a a. 

The isotopic 

distribution of these fragments is governed by the free energy at the density and tempera- 
ture of the emitting system. We have recently discussed the isotope production in terms of 
the Modified Fisher Model |3j. The experimental results exhibit a dependence on an order 
parameter m = N ~ z . This analysis and the often reported observation of isoscaling for 
products of two similar reactions with different neutron to proton ratios, N/Z (3 7 1 make it 
clear that terms in the free energy which are sensitive to the difference in neutron and pro- 
ton concentrations are very important in the fragment formation process. Indeed, isoscaling 
analyses based on the comparison of isotope yields from excited systems of similar temper- 
atures and Z have been employed to obtain information on the symmetry energy and its 



density dependence 



-|8j. The ratio of the isotope yields, R42, between two similar reaction 



systems with different N/A ratios can be expressed by the following isoscaling relation 



R 12 (N, Z) = Cexp(aN + pZ), (1) 

where the isoscaling parameters, a = (/i* — nl)/T and (3 = (/i^ — /ip)/T, represent the 
differences of the neutron ( or proton) chemical potentials between systems 1 and 2, divided 
by the temperature. C is a constant. 

In terms of a modified Fisher model description the experimental yield of an isotope with 
N neutrons and Z protons can be written as [l sj]: 

Y(N, Z) = Y A- T exp{-[G(N } Z) 

-HnN-v P Z]/T}, (2) 

where Y is a constant, G(N,Z) is the nuclear free energy at the time of the fragment forma- 
tion, \x n and \i p are the neutron and proton chemical potentials, and T is the temperature 
of the emitting source. The factor, A~ T , originates from the entropy of the fragment 
Notice that in the ratio of yields employed in an isoscaling analysis the power law term from 
Eq.(E|) will cancel out. 

n a grand canonical treatment the relationship between the isoscaling may be expressed 
as [6|, |7| : 

a(Z) = AC sym A(Z s /A s ) 2 /T, (3) 



where A(Z S /A S ) 2 = (Z s /A s ) 2 1 — (Z s /A s ) 2 2 , C sym is the symmetry free energy and T is the 
temperature. The symmetry free energy is presumed to be density dependent. In a similar 
fashion (3(N) can be expressed as : 

f3(N) = 4C sym A(N s /A s ) 2 /T } (4) 

where A(N S /A S ) 2 = (N s /A s ) 2 1 — (N s / A s ) 2 2 . This suggests that for two systems at tempera- 
ture T in which the symmetry energy is the dominant factor in determining the yield ratios 
of the emitted fragments, the ratio (3(N)/a(Z) can be expressed as 

0(N) (N s /A s )l - (N s /A s )l 
77 a(Z) (Z s /A s )l - (Z s /A s ) 2 ' 1 } 

This ratio will approach -1 when systems of very similar N/Z are considered but in general 
can be quite different from -1. 

It is well known that nucleon forces are isospin invariant, because of this we expect that, 
in the absence of spontaneous symmetry breaking, (/i* — /i^) = (/i* — fi 2 ) and hence 

a = -13. (6) 

In a manner similar to the case for mirror nuclei, at low excitation energies we can expect 
this invariance to be broken by Coulomb energy contributions. However in fragmentation 
reactions occurring near the critical point of the liquid-gas phase change the nuclear symme- 
try should be restored because of the invariance of the nuclear Hamiltonian when m — > —m. 
In a separate paper discussing the analysis of the same data set used in this paper, a clear 
fragment Z (or N) dependence of a{Z) (or 0(N)) has been reported 0, ^From the de- 
tailed comparisons to the dynamical model (AMD) calculation followed by the statistical 
decay code (Gemini), this Z dependence is attributed to the statistical secondary decay 
process of the excited fragments after they are formed at freezeout of the emitting source. 
A significant modification of the isoscaling parameters is also suggested during the cooling 
process. Similar results have been also reported in ref. In that work it is also concluded 
that secondary decay effects play a significant role in determining the observed ratio. Their 
conclusion is based upon the application of an SMM model description to the experimental 
data within a statistical model description of fragmentation. In this paper we suggest that 
there may be an essential relatioship between the scaling parameter and N/Z of the emit- 
ting system, which is related to the restoration of symmetry near the critical point of the 



emitting source and which is sustained in the experimental observables during the cooling 
process and experimentally manifested in the isoscaling parameters. 



I. EXPERIMENTS AND ANALYSIS 

Using high resolution detector telescopes with excellent isotope identification capabili- 
ties, we recently studied a number of heavy ion reactions to determine relative yields for 
production of a wide range of isotopes 0, [9] . The experiment was performed at the K-500 
superconducting cyclotron facility at Texas A&M University. 40 A MeV 64 Zn, 70 Zn and 64 Ni 
beams irradiated 58 Ni, 64 Ni, 112 Sn, 124 Sn, 197 Au, and 232 Th targets. Intermediate mass frag- 
ments (IMFs) were detected by a detector telescope placed at 20 degrees relative to the beam 
direction. The telescope consisted of four Si detectors. Each Si detector had an area of 5cm 
x 5cm. The thicknesses are 129, 300, 1000 and 1000 fim. Using the AE — E technique we 
were typically able to identify 6-8 isotopes for a given Z up to Z=18 with energy thresholds 
of 4-10 A MeV. More details of the analysis are contained in ref. js, 9]. 

Isoscaling analyses were carried out for all possible combinations of these reactions. Eigh- 
teen different reactions are considered here and therefore more than 150 combinations are 
studied. Data for each atomic number were independently fit to extract the isoscaling pa- 
rameter a(Z). /3(A) values were also extracted for each neutron number. For some systems 
the extracted a{Z) parameter shows a steady decrease as Z increases. The /3(A) parameter 
generally showsa much smaller variation with increasing N, and has the opposite sign. A 
clear correlation between them, i.e. a(Z) ~ —/3(A) for the equivalent number of nucleons, 
N = Z, has also been observed as suggested in the introduction(see Eq.([6])). In Fig.l, the 
extracted isoscaling parameters for the case of Z = N = 7 are shown as a typical example. 
Similar correlations are also observed for other selections of Z and N values if Z = N. 

In part a of Figure 1 we have plotted /3(A) vs a(Z). As seen in the figure, the relation 
a(Z) ~ —/3(A) is observed for a(Z) < 0.5 and may deviate slightly at the larger a(Z) 
values. Those larger values are associated with the largest N/A values for the compound 
system. In the bottom part of the figure, these values on the left are compared to predictions 
of Eq.([5]) on the right with the assumptions that Z/A is that of the compound system. We 
note that the experimental values tend to be significantly closer to -1 than the calculated 
values. Except at the low experimental values of a{Z) where the scatter is significant, the 
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FIG. 1: (a) /3(A r ) vs ot(Z) is plotted for N=Z=7. Line indicates the locus f3(N) = —a(Z); (b)ratio 
of f3(N)/a(Z) data to theoretical Eq.© (open squares); (c)/3(N)/a(Z) vs. a(Z) from data (open 
triangles); (d)the analytical prediction Eq.© (open circles) to compare with data. 

experimental values for (3(N)/a(Z) are about 20% lower in the absolute value than the 
model values as indicated by the ratio of these two quantities also plotted in the middle 
part of the figure. In order to see the system dependence of ot(Z) and 0(N) values, these 
values are plotted for separate groups of fissility values in Fig. 2. The fissility is defined 
as X = where Z s and A s are the charges and masses of the source which we assume 
to be the compound nucleus for simplicity. We can define a combined fissility parameter 
between reactions (1) and (2) as EX = X1 + X2. Larger (absolute) values of ce(Z) and /3(iV) 
correspond to large values of the EX parameter. In the figure a(Z) and (3(N) values are 
separately plotted for four different ranges of fissility group for the same data set used in 
Fig.l. We see no systematic correlation with the fissility parameters in the deviation from 
oi(Z) = —(3(N), which might be suggestive of the fact that the Coulomb force is not so 
effective for breaking the invoked invariance of the Nuclear Hamiltonian. It would be very 
interesting to see if Coulomb effects become more important for heavy colliding nuclei such 
as U + U. One should note that a similar result is observed for a(Z) and /3(N) in other 
IMFs, when Z and N are the same. One can also use a and (3 values averaged over a range 
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FIG. 2: /3(N) vs a(Z) with 6 < Z < 13 for different ranges of the fissility parameter. T,X < 60.9 
(top left), 60.9 < Y,X < 72.2 (top right), 72.2 < < 83.4 (bottom left) and 83.4 < Y,X < 94.7 
(bottom right). 

While the trend in FigsfT] and [2] is interesting, it is important to note that when the 
neutron and proton concentrations of the initial excited source are different, two well estab- 
lished trends can act to shift the balance toward symmetric matter and hence to bring the 
absolute values of the observed a(Z) and P{N) parameters closer together. The first is the 



distillation effect in which early emission of particles favors neutron emission over proton 
emission [13I, 



3- 



As a result of this early emission the fragmenting system will tend to 



have a higher symmetry than the initial system. The second is secondary decay of init ially 



rr 



excited fragments |l5|, Il6| which favors a shift toward the evaporation attractor line |17| . 
Thus even if the comparison of primary fragment yields would lead to a significant difference 
in the two isoscaling parameters the subsequent decay can reduce this difference. 
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m-SCALING 



Pursuing the question of phase transitions, we note that we have previously discussed 
some of the present yield data within the Landau free energy description {^J . In this approach 
the ratio of the free energy (per particle) to the temperature is written in terms of an 
expansion: 

F1 2 1 4 1 6 H 

— = —am H — bm H — cm — m— , 7 
T 2 4 6 T y ' 

where m is the order parameter, H is its conjugate variable and a — c are fitting parameters. 
In our case m = (I /A). Notice that the free energy that we have indicated with F includes 
the chemical potential of neutrons and protons i.e. AF(m,T) = [G(N,Z) — jj, n N — \i v Z\ 
(compare to Eq.(j2])). 

We observe that the free energy is even in the exchange of m — > —m reflecting the 
invariance of the nuclear forces when exchanging N and Z. This symmetry is violated by the 
conjugate field H which arises when the source is asymmetric in the chemical composition. 
We stress that correctly m and H are related to each other through the relation m = — 4jr. 

An immediate consequence of the application of the Landau expression of Eq.((7j) in the 
Modified Fisher Model is that it brings a scaling law for m=0 isotopes. Since F(m=0,T)=0, 
for any T, the yield in Eq.(T2]) is given as 

Y(N,Z) = Y A- r , (8) 

for all reactions. In Figj3j yield ratios for m=0 isotopes are separately plotted as a function 
of A for even-even (top) and odd-odd (middle) isotopes for all 18 reactions studied here. In 
order to eliminate the effect of the constant in Eq.(j2J), which are slightly different in each 
reaction system, the yield is normalized to that of the 12 C in each reaction. As seen in the 
figure, the yields from the different reactions are indeed scaled well with A up to A ~ 30 when 
even-even and odd-odd isotopes are plotted separately. One should note that data points 
for a given A represent all 18 different reactions in the figure. The slope difference between 
even-even and odd-odd isotopes can be naturally attributed to the pairing effect. However 
a large pairing effect is expected only at a low temperature, because it is related to the shell 



at a 


high 


11 


Jl2|. 



According model simulations which they have performed the experimentally observed pairing 
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FIG. 3: Yield ratio of m=0 isotopes vs A. The yield is normalized to that of 12 C. Data from all 18 
reaction systems studied in this experiment, (top) even-even isotopes, (middle) odd-odd isotopes 
are plotted, (bottom) pairing corrected yield, Y(N,Z)/exp(<5), are plotted for all m=0 isotopes. 
Lines in each figure are linear fitted ones, r values are 3.3, 2.2 and 2.8 from the top to bottom. 
a p /T =2.2 is used in the bottom. 

effect is attributed to the last chance particle decay of the excited fragments during their 
cooling;. This hypothesis is also supported by our model simulations presented in a separate 



paper 



9|. In order to take into account the pairing effect, data for even-even and odd-odd 



isotopes were simultaneously fitted by the following equation, 

Y(N, Z)=Y A T exp(6/T) 1 



(9) 



8(N,Z) 



(10) 



a p /A l l 2 (odd-odd) 
(even-odd) 
— dp/ A 1 / 2 (even-even). 

and the parameters r and a p /T values was extracted. Using these extracted parameters, 
the experimental yield was divided by the exponent in Eq.()9]) as factor. The results are 
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plotted in the bottom of the figure for all isotopes with m=0. The extracted r value is 2.8 
which is larger than the normal critical exponent 2.3. This difference may reflect either the 
temperature of the emitting source is below the critical temperature or that the secondary 
decay processes modify the value. 

Because of the symmetries of the free energy when we take the ratio between two different 
systems, presumably at the same temperature T and density p, all even order terms in m 
cancel out while the odd terms remain. Those terms depend on the external field H/T. 
Taking the ratio between two systems as in Eq.flSJ) we easily obtain: 

R 12 {m) = Cexp(AH/TmA), (11) 

where AH/T = H\/T —H 2 /T . We can fix the constant C by dividing each experimental yield 
by the 12 C yield following in ref. The goal is to get C — > 1 for reasons that will become 
clear below. Comparing the latter equation with Eq.flT]) we obtain: AH/TmA = aN + (3Z 
i.e. a = —(3 = 4p. As shown in Figfj], for the comparison for isotopes of a given Z with 
isotones having N equal to that Z, this relation appears to be satisfied. The relation is valid 
more in general, and in fact we could write the chemical potentials of neutrons and protons 
as: 



/i n N + fipZ = fiA + HmA, (12) 
from this relation it follows that: 

2H = /i n - n P \ 2/i = n n + n P . (13) 

All these relations show that if m is an order parameter then a = —ft = AH/T. 

The external field is given by the difference of chemical potentials between neutrons and 
protons of the emitting system as expected. ^From Eq. fllll) we can obtain the difference 
between the free energies (or alternatively the external fields) as: 

-ln(R^(m)) = AH/Tm + cQnstant (M) 

Thus a plot of —ln(R)/A versus m should give a linear relation whose slope is given by 
AH/T. This linear relation is demonstrated in Figsfj] and [2] where such a plot is obtained 
for different colliding systems for the isotopes in the selected range of Z. In thateiven range 

nn n n 

a{Z) increases about 50% on average [8j, [9|. As discussed in references [8|, |9j, the observed 
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FIG. 4: Experimental ratios vs m for isotopes with 6 < Z < 13 for (a) tu§^wt^, (b) W^sa ^" , 
(c)S. (d) ^g^ , (e) ^S and ( f ) ^zn"'fn respectively at 40MeV/A The lines are 



, 64 Aff 197 Au 

the results of a linear fit according to Eq.([I 

fragment Z (or N) dependence of the isoscaling parameters is mainly established during 
the statistical cooling of the excited fragments. In fact it has been demonstrated that 
ce(Z) parameter extracted from the primary fragments of the AMD simulations shows no 
significant fragment Z dependence. It should be noted that it is important to normalize the 
distribution ( for instance to 12 C ) as we have done in order that the normalizing constant 
in front of the yield in Eq. fllll) is one. If not this will carry a -j term which might violate 
the scaling. Overall the scaling is satisfied for this set of data as seen in FigJU Compared 
to 'traditional' isoscaling where a fit is performed for each detected charge Z (or each N) 
we see that all the data collapse into one curve. 

We can further elucidate the role of the external field H/T writing the Landau expansion 
and 'shifting' the order parameter by m s which is the position of the minimum of the free 
energy. Such a position depends on the neutron to proton concentration of the source 
Thus 

F 1 1 1 

— = -a{m - m s ) 2 + -b(m - m s ) A + -c(m - m s f. (15) 
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Comparing to Eq.Q we easily obtain 

— = (a + bm 2 + cm 4 )m s + ( 1 — bm H — cm 3 )m 2 + O(mt)---, (16) 

T 2m 2 2 

thus H depends on the source isospin concentration though the parameter a, b, c which are 
terms of the free energy. We stress that these terms refer to the free energy and not to the 
internal symmetry energy. If b and c are of comparable magnitude to parameter a, then 
taking terms of a, Eq.( fl6|) can be further simplified as 

^ = -am 8 + ^a^ + 0(mi)..., (17) 

II. RECONCILIATION OF THE TWO APPROACHES 

Standard isoscaling results have been derived under a general grand canonical ap- 
proach {5-7.]. The Landau approach should be equivalent to it under certain conditions. 
Experimentally the b and c values have not been established because all isotopes identified 
in the present data have m < 0.5 except for nucleons. In the case that b and c are of com- 
parable magnitude to parameter a, which is assumed in the derivation of eq. (|TT|) . we easily 
obtain: 

AH 1 

-jr mA = aAm s (N - Z) - -a(m 2 sl - m 2 s2 )A = aN + (3Z (18) 
which introduces a volume term. Equating similar terms we get: 

a = aAm s — -a{m 2 sl — m 2 s2 ) (19) 

where Am s = m s i — m S 2. It is straightforward to demonstrate the equivalence of the last 
equation to eq.Q derived from the grand canonical approach. In particular we get: 

a + (3 = -a(m 2 sl - m 2 s2 ) (20) 

which shows that the two approaches are equivalent and that m is an order parameter if 
a + (3 = i.e. neglecting 0{m 2 ) terms in the external field. In figure ([5]) we plot a + vs. 
2 2l unfortunately the error bars are rather large but we can see a systematic deviation 



m sl —m 



from zero as expected from eq. (l20l) for large differences in concentration. This indicates that, 
at the level of sensitivity so far acheived with data of this type the presence of higher order 
terms in m is difficult to quantify. Thus, within the error bars, m could be considered 



11 




FIG. 5: ot{Z) + /3(N) vs the difference (solid circles) in concentration for the two reaction systems 
for the case of Z=N=7. The dotted line is the result of a linear fit. 

an order parameter when relatively neutron (or proton) rich sources are considered. In 
particular, phase transitions in finite systems could be studied using the same language of 
macroscopic systems i.e., 'turning on and off' an external field jsj. 



SUMMARY 



In conclusion, in this paper we have discussed scaling of ratios of yields from different 
colliding systems under similar physical control parameters, i.e. density and temperature. 
A careful and precise determination of isotopic yields is needed in order to see the features 
of the system near the phase transition. There is an order parameter, m, given by the dif- 
ference in neutron and proton concentrations of the detected fragments which leads to an 
expected isoscaling relation, a direct consequence of the restored symmetry of the nuclear 
Hamiltonian when exchanging neutrons with protons. The data suggest that the Coulomb 
field may not significantly violate such a symmetry. The existence of m-scaling might be a 
signature for near criticality of the fragmenting system. Other properties of the 'rich' nu- 
clear Hamiltonian, such as pairing, appear to result in small violations of the scaling. This 
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is an interesting physical aspect which deserves further and deep investigation both theo- 
retically and experimentally Also, it would be interesting to search for m-scaling violations 
in heavily charged colliding systems such as U+U. The absence of a violation in these cases 
would suggest that densities and deformations of the fragments are such that the effect of 
Coulomb is significantly reduced. Studies of the other extreme case of very exotic colliding 
systems would be also be valuable to probe the effects of high 'external' field on the phase 
transition. The atomic nucleus constitutes a formidable laboratory to test our knowledge 
and understanding of phase transitions in a finite system and offers a unique possibility for 
different quantum aspects similar to other bosons and fermion mixtures. 

A major consideration in the interpretation of the results presented in this paper is the 
effect of the secondary decay process. In the experiments excited fragments cool down to the 
ground state before they are detected. The reconstruction of the primary fragments from the 
experimentally observed IMFs and associated particles is not straightforward, since multiple 
excited primary fragments may be simultaneously produced in multifragmentation reactions 
making the unambiguous identification of the primary fragment distribution difficult. Indeed 
a major goal of the experiments from which the present isoscaling data are taken was to 
employ fragment-particle correlation measurements to reconstruct the primary fragment 
distribution. The correlation data are still being analyzed jiol ]. 
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